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Abstract: We consider the most general supersymmetric model with minimal particle 
content and an additional discrete Z% symmetry (instead of R-parity) , which allows lepton 
number violating terms and results in non-zero Majorana neutrino masses. We investigate 
whether the currently measured values for lepton masses and mixing can be reproduced. 
We set up a framework in which Lagrangian parameters can be initialised without recourse 
to assumptions concerning trilinear or bilinear superpotential terms, CP-conservation or 
intergenerational mixing and analyse in detail the one loop corrections to the neutrino 
masses. We present scenarios in which the experimental data are reproduced and show 
the effect varying lepton number violating couplings has on the predicted atmospheric and 
solar mass 2 differences. We find that with bilinear lepton number violating couplings in the 
superpotential of the order 1 MeV the atmospheric mass scale can be reproduced. Certain 
trilinear superpotential couplings, usually, of the order of the electron Yukawa coupling can 
give rise to either atmospheric or solar mass scales and bilinear supersymmetry breaking 
terms of the order 0.1 GeV 2 can set the solar mass scale. Further details of our calculation, 
Lagrangian, Feynman rules and relevant generic loop diagrams, are presented in three 
Appendices. 
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1. Introduction 

Once one allows for R-parity [1] violation in the Minimal Supersymmetric Standard Model 
(MSSM) there is an embarrassingly large class of possible models. Building on the sem- 
inal paper of Ibanez and Ross [2], it has been shown that there are three preferred Z^ 
symmetries which can be imposed when constructing the MSSM with minimal content of 
particle fields [3]: one is the standard R-parity under which the Standard Model particles 
are R-parity even while their superpartners are R-parity odd, the second is a unique, Z3 
symmetry which results in the MSSM with lepton number violation; denoted here as ( 
JZ-MSSM ) and the third is a Zq symmetry refered to as proton hexality. The former guar- 
antees a stable lightest supersymmetric particle and thus missing energy at colliders, the 
Z3 and Zq lead to proton stability 1 together with neutrino flavour changing phenomena 
and masses. In this paper we want to investigate in detail neutrino masses and mixings in 
the 1/i-MSSM motivated from the current observations of neutrino flavour metamorphosis 
and proton stability [4]. 

In the 1/i-MSSM a single neutrino mass arises at tree level due to the mixing between 
neutrinos, gauginos and higgsinos [5-8]. This tree level mass is proportional to the bilinear 
lepton number violating superpotential parameter, /ij, squared, which is assumed to be of 
order 2 of MeV, and is suppressed by the "TeV" supersymmetry breaking gaugino masses, 
resulting in a low energy see-saw mechanism with light neutrino and heavy neutralino 
masses. The other two neutrino masses arise from quantum loop corrections made up from 
lepton number violating superpotential or supersymmetry breaking vertices. We shall refer 
to these neutrinos with the term "massless neutrinos" . 

Calculations for neutrino masses in the JZ-MSSM have been addressed many times in 
the literature. The tree level set-up of the model was first given in [5], and details worked 
out later in [6-8]. Calculations of the one-loop neutrino masses including only the bilinear 
superpotential term are given in [11—15]. Corrections involving the trilinear superpotential 
Yukawa couplings A, A' considered mostly in the mass insertion approximation [17-20] and 
under the assumption of CP-conservation and flavour diagonal soft SUSY breaking terms. 
Renormalization group induced corrections to neutrino masses have been studied in [21-23]. 

1 Baryon number violating operators of the form QQQL or UUDE are not allowed in I/-MSSM . 
2 This can be naturally accommodated by employing an R-symmetry [9,10,23]. 
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There is of course a vast number of articles using these calculations, or simplified versions 
of them, in order to describe the solar and atmospheric neutrino puzzles [24, 26-34] . 

In this work, we calculate the complete set of the one-loop corrections to the massless 
neutrinos without resorting to approximations about CP-conservation or bilinear superpo- 
tential operator dominance. The outline of our paper is as follows : in section ^ we show 
how to define the Lagrangian parameters in the fermion sector of the theory, by starting 
out with physical input parameters like the lepton masses and mixing angles. In section || 
we describe our renormalization procedure and present analytical results for the one-loop 
corrections together with approximate formulas (if necessary) for individual diagrams. We 
compare with the current literature. In section |I] we present numerical results for the size 
of the input parameters when these account for the neutrino experimental data, and in 
section ||] draw our conclusions. Finally, in three Appendices, we set out our notation, in 
Appendix |A|, the Lagrangian and the mass terms, and present the relevant Feynman rules 
in Appendix [B|. Furthermore, in Appendix ^ we present a pedagogical brief introduction to 
the Weyl spinor calculation that are employed throughout this paper and present general 
one-loop self energy corrections that are employed in this paper and can be used elsewhere. 

The Fortran-code for calculating the neutrino masses used in this article has been 
made publicly available 3 . It can be used in adding additional constraints when other 
L'-MSSM processes are studied. 

2. Fermion masses and mixings in ^-MSSM 

Fermion masses and couplings of the general L'-MSSM are defined by the superpotential 
of the model, vacuum expectation values (vevs) of the neutral scalar fields and the soft 
supersymmetry breaking gaugino masses. The most general superpotential takes the form 

Wl/_MSSM = e ab (l^^afik-C-a^Ek + ^'ajk^aQ^ Dk,x + (Yu) jkQ°j' X ' H\lJk,x ~ 

(2.1) 

where Qf x , Df, Uf, Cf, E{, Hf, H% are the chiral superfield particle content, i = 1, 2, 3 is 
a generation index, x = 1,2,3 and a = 1, 2 are SU(3) and SU (2) gauge indices, respectively. 
The simple form of ( |2.1| ) results when combining the chiral doublet superfields with common 
hypercharge Y = — \ into £^=0... 3 = (Hf, Lf =123 ). fi a is the generalized dimensionful 
//-parameter, with hq and = 1, ...3 the lepton number conserving and violating parts 
respectively, and \ a (3k, ^' a jk> (Yu)ij are Yukawa matrices with e a b being the totally anti- 
symmetric tensor ei2 = +1. 

Physical masses of the fermion fields depend on appropriate A, A',// and Yjj couplings 
multiplied by the vevs of the neutral scalar fields. As it has been shown in [35,36], by unitary 
rotation in the 4-dimensional space of the neutral scalar components of C a , it is possible 
to set three of the four vacuum expectation values of the C a fields to zero, leaving two real 



3l 
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non-zero vevs in the neutral scalar sector and, simultaneously, significantly simplifying its 
structure. It is convenient to apply such a transformation not just to scalars, but to the 
whole chiral super field, 




The tildes and primes on the fields are then dropped. 

In a standard way, both isospin components of Q superfield and of U, D superfields 
can each be redefined by a unitary rotation in the flavour space. As such, it is possible to 
diagonalise the Yukawa couplings (Yd), (Yrj) (note that (Yd) = A^- in the basis with 
two non- vanishing scalar vevs) and absorb the rotation matrices in field redefinitions such 
that they do not appear explicitly in the Lagrangian, apart from a specific combination of 
rotation matrices which appear in the gauge and Higgs charged currents which is identified 
as the CKM matrix. In this basis, it is clear how to initialise the Lagrangian parameters, as 
the diagonal values are then proportional to the measured values for the up- and down-type 
quarks. For more details concerning this point the reader should consult Appendix [A]. 

In the lepton sector, however, the same approach cannot be adopted for two reasons. 
Firstly, even with a diagonal Yukawa matrix, the charged lepton masses are given by three 
eigenvalues of the larger (5 x 5) mass matrix which includes mixing between the charged 
fermionic components of the £, E and the charged gauginos and higgsinos. Thus, the 
diagonal entries in the Yukawa matrix would not correspond exactly to the masses of the 
physical mass eigenstates which describe the charged leptons. Secondly, the £-basis has 
already been fixed by the property that three neutral scalar vevs should be zero, so we are 
not free to absorb a rotation matrix 4 . 

Still, there is some freedom remaining due to the fact that the .E-base has not, as yet, 
been fixed. Flavour rotation in the E-space can be used to remove some of the unphysical 
degrees of freedom in Yl = Xoij coupling 5 . As every general, complex matrix, Aoij can 
be uniquely decomposed (polar decomposition theorem [38]) into a product of positive 
semi-definite Hermitian matrix Aq and unitary matrix Ve '■ 



basis and to diagonal Yukawa couplings [37]; it is possible to use the freedom in the 3-dimensional lepton 
space, which we used in Ref. [35] to diagonalise the sneutrino masses, in order to make the lepton Yukawa 
couplings diagonal. But then one will have a 10 x 10 mass matrix for the neutral scalars because this 3x3 
rotation is, in general, complex (unitary). We want to avoid this complication by all means. 

If the decomposition is unique, then all unphysical degrees of freedom will be removed, because then 
the full U(3) rotation is absorbed into E and every rotation in the £ a -space will "destroy" some of the 
properties we want to keep. 




4 This is not entirely true: it is actually possible to perform a rotation into the vanishing sneutrino vev 
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Ve can be then adsorbed in the chiral superfield E redefinition and the 'hat' over A is 
dropped. 

After all the transformations described above, we arrived to the form of the superpo- 
tential (pj] ) where (Yd)^ = Aq^- and (Yrj)jj are flavour-diagonal and (Yl)^- = Xoij is 
hermitian. Other coupling constants are free and, in general, complex parameters. 

2.1 Block Diagonalising 

In the following sections we outline the procedure by which the parameters of the general 
hermitian matrix (Yl),„- can be initialised such that the correct values are obtained for 
the charged lepton masses and the MNS mixing matrix [39]. In order to do that, it will 
be convenient to diagonalise the neutralino-neutrino and chargino-charged lepton mass 
matrices in two stages. First an approximate, unitary or biunitary transformation will 
result in matrices in block diagonal form; the standard model and supersymmetric fermion 
masses being split into separate blocks. Then, a second transformation will diagonalize the 
blocks. 

The block diagonalisation can be performed for any complex matrix. Every general 
matrix M can be diagonalised by two unitary matrices V, U : 



V^MU = M = dmg{m 1 ,m 2 ,--- ,m n ) = 1 , (2.5) 




where mf are eigenvalues of MM^ and Mi , M 2 are two diagonal sub-matrices of a chosen 
size. Hence, one can always rewrite M in the form 

M = VMU ] = VA ] AMB ] BU ] , (2.6) 

where A and B are some unitary matrices of the form 



A = \ n 1 ; , B = „ „ , (2-7) 





with sub-matrices A\^, £?i,2 which are also unitary. Thus we can write 

M = QM B P^ , 

where 



M B = AMBi = [ AlMlB i ° t ] , (2.9) 
1 A 2 M 2 B\ 1 



is block diagonal in form and Q = VA\ P = UB^ . Of course, Mb is not uniquely defined. 

Block diagonalisation is particularly useful in case of hierarchical matrices, when one 
can find analytical approximate formulae for P, Q matrices in eq. . Consider for exam- 
ple a hermitian matrix (other cases can be considered analogously) of the form: 

M=( m f mB ), (2.10) 



B m C 
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where itia = m A , mc = m c an d ll 771 ^!! 3> H m s||) ll m c||- In such case, the approximately 
(up to the terms Odlm^cHVHm^H 2 )) unitary matrix U, 

transforms M into approximately block-diagonal form: 

V ] MU _ ( m A+{m A T 1 m B m ] B + m B m ] B (m A y 1 (m A y 1 m B m c \ | g f \\m B ,c\\ 3 
\ rr^ c m) B (m a)~ X mc — m B (mA)~ x m B J \ \\mA\\ 2 

° ^+o( llmB '€) ■ (2.12) 



mc — m B (m A ) 1 m B J V \\ m A 



We kept explicitly the O ( i? 5 |i ) term in (22) element of block-diagonalised form of 



\\ m A\\ 

M as in many models mc = and in this case it will be the only term which survives 
(the see-saw mechanism [40,41] for neutrino masses being the most famous example of this 
hierarchical structure). 

As a next step one needs to find matrices that diagonalise the sub-blocks in eq. (|2.12j ). 
We employ this method in the next sections where we present explicit perturbative (l s t 



order) results for analogs of the matrices Q and P in Eq. (2.8), in both neutral and charged 
fermion masses in the JZ-MSSM . In passing that, although we use only the approximate 
analytical expressions for the see-saw type expansions above, in our numerical predictions 
we perform exact block diagonalization, iteratively finding the correct, and strictly unitary, 
matrices P, Q of eq. (E 



2.2 Fermion masses and mixing 

In the following section we shall present the tree level phenomena of the fermion sector 
in the 1/i-MSSM . We consider in turn, the neutral and charged fermion sectors and the 
patterns of the mass matrices. We consider the tree level eigenvalues, particularly for the 



neutral sector and the approximate block diagonalisation of the matrices. In section ^3 
we use the approximate block diagonalisation to consider the way in which the MNS matrix 
appears in this model, as it is now a sub-block of a larger unitary matrix, the MNS matrix 
itself is not generally unitary. This analysis is then used to ensure the correct low energy 
parameters are reproduced, despite mixing between the leptons and heavy fermion fields. 

In section [3] we consider the effect of radiative corrections at the order of one-loop. After 
setting out the renormalisation framework, we present in turn various loop diagrams and 
highlight the important contributions. The full numerical analysis has been completed, 
however approximate expression are presented for each contribution which demonstrate 
from where the important effects arise. We will consider the case where the tree-level effect 
dominates and gives rise to the larger, atmospheric mass squared difference, in which case 
the solar mass squared difference is generated by the loop effects. We also show that it is 
possible for loop effects to be greater than the tree level affects, in which case both mass 
squared differences are generated at the level of one-loop. 



-6- 



2.2.1 Neutral fermion sector 

In the lepton number violating extension of the minimal supersymmetric standard model 
( ^-MSSM ) the neutrinos (^123)' neu t ra l higgsinos (v L0 and h^and neutral colourless 
gauginos (W and B) mix. To transform the fields into the mass basis, the 7x7 neutralino 
mass matrix must be diagonalised. In the interaction basis, 



CD--[-iB, -iW , h° 2 , y La )M N 



( -iB \ 
-iW 
hP 2 

\ U L(3 J 



where the full 7x7 mass matrix reads 



+ H.c 



(2.13) 



a a I Mn4 x 4 <ijV4x3 \ 

Mn = \ n 

\ "JV3X4 U 3x3 J 

and the sub-blocks are, in the basis (—iB, —iW , h®, v La = h®, vi) [35] 



(2.14) 



M N 



4x4 





Mi 





gvu 
2 


gv d 
2 


\ 







M 2 


2 


92V d 
2 






gvu 
2 


g2v u 

2 





-Mo 




V 


2 


92 v d 
2 


—fJ>o 





/ 



(2.15) 



and 



dN4 



x3 



V 






-Ml 















"M3 





(2.16) 



There is no quantum number to differentiate between neutralinos and neutrinos, the states 
of definite mass do not have definite lepton number and, as such, there is no reason to think 
of neutrinos and neutralinos separately. However, for realistic values of parameters, four 
of the mass eigenstates are heavy and three are very light, so it is convenient to refer to 
them as to neutralinos and neutrinos, respectively. In addition to this, it can be seen that 
the mixing is sufficiently small that these three light neutral states are the states which 
dominantly appear in the decay of the W boson to charged leptons, differentiating between 
the eigenstates we refer to as neutrinos from those we refer to as neutralinos. 

The matrix, Z^;, which rotates the fields in ( 2.13j ) from the interaction basis to the 
mass eigenstate basis is given by 



/ 



-iB \ 
-iW° 
h° 2 



Zn 




(2.17) 



where 



\ "La J 

7 are seven neutral two component spinors 
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The matrix A4n, as it has been split in eq.( 2.14| ), contains block diagonal terms that 
conserve lepton number and off diagonal blocks which violate lepton number. The latter 
are expected to be very small, as they are strongly constrained by the bounds on neutrino 
masses or other lepton number violating processes. Thus, one can use the block diagonal- 
ization procedure of section and, neglecting terms of the order j£- and assume Zn to 
be of the form 



'■<N 



d N M V 




The first matrix on the RHS of eq. ( 2,18] ) which is the analog of the matrix Qt in 
diagonalises the neutrino-neutralino mass matrix: 



1 

-d T N M~ l 



M 



1 



N 



M N 



4x4 




(2.18) 
) block 

(2.19) 



1 J v d * M * 1 / V 

where the "TeV" see-saw suppressed effective 3x3 neutrino mass matrix is given by [6-8] 



«// JTM-ij v d(Mig% + M 2 g 2 
" = ~ d » M » dN = ABet[M N ] 



' (4 MlA*2 MlM3\ 
MlM2 l4 /^M3 



(2.20) 



Physical neutralino masses and mixing matrix Zn can be found in a standard manner by 
numerical diagonalization of the matrix Mpf. Diagonalization on mt can be easily done 
analytically, leading to two massless and one massive neutrino, with its mass given by: 



in 



tree 



v 2 d {M ig l + M 2 g 2 



4Det[M J 



N 



and the mixing matrix Z v is 



Z„ 



VM 2 +M 2 



ImiIIwI 



(M 2 + IH 2 + I M3 1 



Iwl 



(2.21) 



^HW^WlW vVil 2 +IH 2 +lwi 2 



/WlwP+lwP ImiIVImiP + I^PVImiP + I^P + ImsI 2 /4VM 2 + M 2 + M 2 



V 







mi 



IHVM 2 +IH 2 



-^2x2 








-1 



M3M\/M 2 +IH 2 +|Af3| 2 MiVM 2 + IM2| 2 + |w| 2 / 



(2.22) 



where X 2 ^2 is an SU(2) rotation. At tree level, the five massive eigenstates are unam- 
biguously defined by diagonalising the mass matrix. The two massless eigenstates, due 
to the fact that they are degenerate in mass, are not fully defined. The eigenstates are 
chosen to be orthogonal, but it is still possible to perform a rotation on the eigenstates. 
As such, statements about the lightest neutrinos, v\p, are basis dependent. Because of 
this, the one loop contributions to M.Npq are also basis dependent. By choosing a different 
linear superposition of the tree level eigenstates, the one-loop contributions to the 2x2 
sub-lock M. N (5 5 6)(5,6) referring to the massless neutrinos would be redistributed between 
themselves. This freedom of basis choice is only present at tree level and is not physical. 
Thus, we start from X 2 x2 = 12x2 and after calculating the radiative corrections to the 
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neutralino-neutrino mass matrix we adjust X 2x2 such that the off-diagonal one- loop con- 
tribution 5Mn 56 is approximately zero (this can be done iteratively) . As such the effect of 
rediagonalising the neutrino sector after loop corrections are added is small. As we discuss 



in section 2.3, choosing the basis in this manner helps also to define the lepton Yukawa 
couplings in terms of measured quantities like lepton masses and the Umns mixing matrix. 

The result that two of the neutrino masses vanish at the tree level is not the effect of 
the approximations made [6-8] . The explicit calculation of the secular equation for the full 
neutralino-neutrino mass matrix A4 at , results in 



det(A^Ar - A) 



-A 2 [Adet(MTv-A) 
(/x 2 + f4 + (4) (A (Mi 



A)(M 2 -A) + ^(Mi 



A) + ^(M 2 



(2.23) 
A)' 



Hence, M n always has at least two zero modes. This can be seen directly, by noting that 
the final three columns of the 7x7 mass matrix are proportional to each other. 

Finally, the physical eigenstates of neutralinos and neutrinos are approximately given 
by, respectively: 

/ -iB \ 



-iW 



'LO 



+ M N l d N v Ll 



(2.24) 



and 



f4 



\4. 



-4M+T 1 



/ -iB \ 
-iW 

h°2 



\ 



(2.25) 



J 



with v L0 = h*l is the down type Higgsino. For a quick view of definitions see also Ap- 
pendix [A]. 

2.2.2 Charged fermion sector 

In a similar fashion to the neutral sector, charged leptons, gauginos and higgsinos mix. 
The full 5x5 chargino mass matrix in the basis of Ref. [35] is given by 



M 



c 



M C2x2 

dc3x2 mc3x3 



with the lepton number conserving sub-locks 

Mc 2x2 : 



M 2 3^ 



92Vd 
V2 



Mo 



Vd , 

m Cij - ^ A 0ij , 



(2.26) 



(2.27) 



and the lepton number violating being 

dc3x2 




(2.28) 
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The rotation matrices which transform between interaction eigenstates and mass eigen- 
states are given by 



(2.29) 









(4\ 




h+ 






V 


e Ri J 









z* 



and, as such, the mass matrix is diagonalised 

M c = Z ] _M C Z+ , 



(2.30) 



(2.31) 



where the 'hat' denotes that the matrix is diagonal. 

The matrices Z + and Z_ can be determined by the requirement that they should 
diagonalise the Hermitian matrices A4 c A4c and A4c-Mq, respectively. The off-diagonal 
blocks in the latter two combinations are small comparing to the diagonal ones, so one can 
again use block-diagonalising approximation of section |2.1| . Keeping just the leading terms 
in 1/Mc expansion, one obtains 



1 



d c M7, 1 



-M c ~% 





(2.32) 



Substitution of (2.32) in ( p. 31 ) results in the physical effective mass matrix 



M 



c 



o 



d C 
Ml 



Zl_mcZ l+ 



(n ( d c m c \ 
U \M C ) 

n ( d 2 c m c \ 



zIm g z^ 









Z _mcZiA 



(2.33) 

Then the matrices Z + , Z- can be again determined as diagonalising matrices for the 
MgMc, McMg products, with the additional requirement that physical fermion masses 
are real and positive. Matrix mc in our basis is hermitian and as such Zi+ = Z\- = Zi. 
Furthermore, physical eigenstates of fermion fields are given by 



h+ 



Z, 



e R2 



(2.34) 



- 10 - 



and 



-iW~ 



-LO 



+ Mc 1T d T c 



I K 3 \ 



-d* c Mc U 



-iW 



e L2 
\ e L3 

(e 



+ 



Li 
VL3 



(2.35) 



For a quick view of the full charged fermion mass matrix see Appendix |A|. 
2.3 Constructing the MNS matrix 

The lepton mixing matrices appear in the charged current gauge boson vertex. Whereas 
in the lepton number conserving case the C/mns matrix is a 3 x 3 matrix describing the 
mixing of three charged leptons into three neutral leptons, the R-parity violating case has 
the mixing of five charged fermions into seven neutral fermions, of which the Umns is a 
3x3 sub-matrix, only approximately unitary. Thus, 



92 



C D ^= W+v' Li a»e' Li + H.c = ^= W+ R° p <t» (U MNS [ 



pq K q + H.C 



where primes refer to interaction eigenstates, and the MNS matrix, 

/ d c d N \ 
\McM N J ' 



C/mns = Ztzf + O ( J^-] 



(2.36) 



(2.37) 



is defined in terms of the mixing matrices introduced in ( 2.22 ) and below fl2.33j ). 

As the first term in Eq. ( p. 37 ) is unitary, unitarity violation in £/mns is at most of 



the order of 



2 M SI 



M 



tan^ [3 ~ 10 12 tan 2 (3, which is well below sensitivity of 



current (or planned) experiments determining the MNS matrix. 
2.4 Input parameters 

The parameters characterising the light charged fermions are already very well known; 
masses are measured with very good accuracy. In contrast to this, the neutrino sector is 
not, as yet, known with the same precision. There is however, information about the mass 
square difference between neutrinos and the mixing between different interaction states, 
and upcoming experiments should improve our knowledge of neutrino parameters in the 
near future. Furthermore, super symmetric fermions have not yet been discovered, and their 
masses and couplings (those which are not determined by supersymmetric structure of the 
model) are entirely unknown. In the JZ-MSSM both sectors mix, and thus the question 
of effective and convenient parameterisation arises. In this section, we will consider the 
parameters in the Lagrangian which effect the tree level masses and mixing. Later, we 
discuss parameters which affect the neutral fermions at the order of one loop. 

As the SUSY sector has not been measured directly, it is convenient to take as an input 
the following set of Lagrangian parameters: Mi, M2, tan/3, (j, = (j,q. With /ij of the order 
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of MeV, corrections to the supersymmetric sector from the light fermion sector are see-saw 
suppressed and negligible. Chargino and neutralino masses and lepton-number conserving 
couplings are thus to a very good accuracy determined by the above four parameters. 
Reconstructing their values from the actual experimental measurements has already been 
discussed in the literature [42]. 

In a next step, neutrino masses can be parameterised at tree level by setting the lepton- 
number violating parameters fii, i = 1,2,3. In the future, when the neutrino mass matrix 
is known to better accuracy, it could become more convenient to reconstruct //j from the 
experimental data - for that, the knowledge of radiative corrections to the neutrino masses 
would be vital. 

To initialize the Lagrangian parameters in the light charged lepton sector, one needs 
to input the lepton masses, m e , m„, m T and the mixing matrix C/mns- The lepton rotation 
matrix Z\ can be then calculated from ( p. 37 ) 



Zi f« Z* v U* MNS . (2.38) 



As we have seen from ( p. 22 ), the neutrino mixing matrix, Z v is defined at tree level up to a 



U (2) rotation for a given set of [/,{. The same matrix is then defined completely at one-loop 
where all the neutrinos are no longer degenerate. Thus, a complete definition of Z[ requires 
a one-loop corrected neutrino mixing matrix. Then, the light charged fermion mass matrix 
mc in ( 2.27| ), which is hermitian and proportional to the Yukawa matrix Aoi? = ^mcij 



is given by: 

m c = Zi diag(m e , m M , m T ) z\ . (2.39) 
Eq. (2.39) holds under the assumption that one-loop corrections to Z[ are small. Otherwise, 



one needs to find mc iteratively, such that physical (i.e. loop corrected) Z v and Z[ produce 
the correct experimentally measured Umns matrix of eq. Q2.38p . 

As we have repeatedly mentioned so far, it is important to notice that the matrix mc 
is not diagonal. 



3. One-loop neutrino masses in ^,-MSSM 
3.1 Renormalization issues 

As we have already seen in the previous section, the presence of the bilinear lepton number 
violating mass term in the superpotential, fii, triggers the mixing between neutrinos and 
neutralinos. Diagonalization of the full 7x7 neutralino mass matrix generates four heavy 
"neutralino" masses and one "neutrino" mass at tree level. Furthermore, the two remaining 
neutrinos become massive at the one loop level due to the presence of other lepton-number 
violating couplings and masses. 

Physical neutralino masses are defined as poles of the inverse propagator. The appro- 
priate formula can be derived from flC-6|) and the definition for the one-particle irreducible 
(1PI) self-energy functions, 
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where the momentum q 11 flows from left to right, and are the Pauli matrices . The 
requirement that the determinant of the inverse propagator is zero, leads to the expression 
for the physical neutralino mass matrix 

mP Npq = r^N^iV-R^pq + [^^N pq (™>Np) ~ m Np pq( m Np)] , (3-3) 

where /jr is the renormalization scale and S^'^ the 1PI contributions to the effective action 
defined ( p|j3.2| ). m^p are the diagonal tree level neutrino masses (they are zero for the 
two massless neutrinos). Our renormalization analysis is similar to the one in Ref. [11]. We 
have also studied an on-shell renormalization analogous to the one in [44]. In this scheme, 
the physical mass formula is similar to (|3.3| ). 

Some additional remarks regarding the present are in order : 

a) The two one-loop induced neutrino masses are perfectly defined at one loop through 
Eq. (|3.3[) . We have proven both analytically and numerically that these masses are 
finite and numerically that they are gauge independent at one-loop order. This result 
remains valid also when one takes into account mixing between them. 

b) The one- loop formula for neutralino masses ( |3.3| ), receives, in addition to diago- 
nal corrections, off diagonal ones, Tn^ e q — > Trij^. Physical neutralino and neu- 
trino masses are then obtained from the diagonalization of Ti^f as rri^ le = (1 + 
5Z^ N )m^ e (1 + 5Zjy). Then the corrected mixing matrix Zn has to be replaced 
by Zn — ► (1 + 5Zn)Zn everywhere in our expressions for the self energies. How- 
ever, the corrections to Zn matrix are of the order 5Z^ pq ~ SmN pq /(mN P p — mNqq) 
and are small if the tree level masses are not degenerate. In our case this happens 
only for the two massless neutrinos, so we include off-diagonal corrections, shifting 
Z v — ► {\+8Z v )Z u where 5Z y has only the upper 2x2 block non-trivial 7 . As discussed 
previously in section [2|, this is actually necessary to fix the neutrino basis. The re- 
sulting corrections to lighter neutrino masses are formally two- loop, but numerically 
important and have to be taken into account. One should note that, as mentioned 
in the previous point, one- loop corrections to the light neutrino mass sub- matrix are 
finite - going beyond the approximation described above would require performing 
formal renormalization on the neutral fermion mass matrix. Finally, similar consid- 
erations apply to the case of the charged fermion mixing matrix Z\ in ( p. 38 ). 

6 We use Weyl spinor notation in our calculation. The corresponding formulae for Weyl-propagators and 
vertices are defined in Appendix ^ and in [45] . 

7 Possible exception is the case when /ii parameters are very small, so that one-loop corrections to the 
neutrino masses are of the order of tree level neutrino mass or bigger. In this case one needs to rediagonalise 
the full 3x3 neutrino mass matrix. This is done numerically in section ^ when presenting our results for 
Mi = 0. 
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c) As we have already mentioned, in our neutral scalar basis [35] the sneutrino vevs are 
zero at tree level. Non-zero sneutrino vevs will appear in general at one-loop. As a 



result, the neutrino tree level mass in Eq. (2.21) should be corrected. However, loop 



induced vev contributions do not arise for the massless neutrinos-they are generated 
outside the 2x2 light neutrino mass matrix-which is the case we are interested in. 



d) We choose = Mz as renormalization scale in ( p.3j ) were we input the DR param- 
eters for m^ e (^iji) at tree level. These parameters are taken after diagonalising the 
full neutralino mass matrix in Eq. ( |2.14| ). 

e) The infinities which arise in the calculation of the one loop corrections, must be 
absorbed in parameters of the tree-level Lagrangian. It is possible to check that there 
are no infinities which must be absorbed where the mass matrix contains zero entries. 
The divergent parts of the integrals do not depend on the masses of the particles in 
the loop integral and as such the infinities only arise when a diagram exists in the 
interaction picture with only a mass insertion on the fermion in the loop. That is, the 
symmetry which prevents a term existing in the classical Lagrangian also causes the 
divergent part to cancel in the mass basis. This guarantees that it is always possible 
to absorb the infinite part of the integral in the bare parameters of the classical 
Lagrangian. 

Having considered the above points we find that for the calculation of the one-loop 
corrections to the eigenstates which are massless at tree level, it is sufficient to consider 
corrections to the bilinear terms purely between these eigenstates. We find that it is 
possible to neglect the one-loop effects which correct other entries of the neutral fermion 
mass matrix, describing neutralino masses of neutralino-neutrino mixing. 

3.2 One loop contributions to the massless neutrino eigenstates 

The mixing of neutrino, neutral Higgsino and neutral gaugino interaction eigenstates has 
been shown to result in two mass eigenstates with zero mass. It is important to note the 
composition of the massless eigenstates; they consist solely of neutrino interaction states, 
not containing any contribution from the fermionic components of the gauge supermulti- 
plets or the Higgs supermultiplets. This can be stated, entirely equivalently as the rotation 
matrix in Eq. ( p. 18 ) becomes [6] 



Z N{l-+4}{5~>6} = • (3-4) 

Radiative corrections at one-loop will affect all three of the light mass eigenstates ('neu- 
trinos') and will lift the degeneracy between the massless eigenstates. The possibility that 
the hierarchy of mass differences in the neutrino sector can be explained in the L'-MSSM is 
considered. If the 'atmospheric' mass difference were to result from the tree level splitting 
and the 'solar' mass difference originated from loop effects, the distinct hierarchy could be 
accommodated within the model. If the solar mass difference is to originate purely from 
loop corrections to massless eigenstates, we must find loop corrections from diagrams with 
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external legs comprised purely of neutrino interaction states. A small caveat is required to 
compare with the literature. In a general basis where the sneutrino vacuum expectation 
values are not zero, the massless neutrinos are comprised of interaction state neutrinos 
and the interaction state Higgsino which carries the same quantum numbers. In the 'mass 
insertion'-type diagrams, this means that only diagrams without mass insertion or with a 
mass insertion which changes the original neutrino external leg to the down-type Higgsino 
can contribute to the solar mass [15], if the assumption that the solar mass arises purely 
from loop corrections to eigenstates which were massless at tree level. 

The one- loop, one-particle irreducible self energies needed in (|3.3| ) are calculated in 
Appendix [C], see ( |C- 1 4 C-To] \C- 1 8 JC- 1 9|) . Results are presented for general vertices and for a 



general gauge. One then has to just replace these vertices with the appropriate Feynman 
rules of Appendix ^ in order to obtain T, D,L . Since this, rather trivial replacement, leads to 
rather lengthy formulae for the self energies, we refrain for presenting the full expressions 
here. Instead we examine in detail the dominant contributions to the massless neutrinos, 
i.e., contributions to T, D . Of course, our numerical analysis exploits the full corrections. 
From the expressions ( P-14 ,C-15), it can be seen that these corrections are propor- 



tional to the mass of the fermion in the loop. As such, the diagrams that give a large 
contribution are the diagrams with sufficiently heavy fermions compared to any suppres- 
sion from the vertices. In addition, standard model neutral fermion masses arise entirely 
due to Supersymmetry breaking in the Fj-MSSM so corrections are expected to be large 
for individual diagrams or a certain amount of fine tuning is required for large SUSY soft 
breaking masses. 

In the next section of this chapter we analyze analytically all the possible contributions 
to for the massless neutrinos, isolating the dominant ones. For simplicity, we shall 
confine ourselves only to the diagonal parts of , although our numerics account also for 
the off diagonal effects in the massless neutrino sub-block. 

3.2.1 Neutral fermion - neutral scalar contribution 

Diagrammatically this contribution reads as: 

✓ / 
/ / 





,.o ,,0 ,,0 



This can be easily calculated by using the formula ( p-14 ) and with the Feynman rules 
read from Appendix [B,l| . The result for the full contribution to the massless neutrinos, 
p = q = {5,6}, is : 



J Npp 



r=l s=l i,j=l 



m K o 
(4^)2 
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2c w 



ZN(4+i)pZNlr 



2s 



-ZN2rZN{A+i) 



^N(i+j)p^Nlr 



2sw 



ZN2rZ]S[(4 + j)p 



Z R(2+i)s Z R{2+j) s Bo(m 2 K o,mj { o,rn 2 K o) - Z A (2+i) s Z A(2+j)sBo(ml,o, m A o, m 2 o) , (3.5) 

where H® 5 and A® 5 are the CP-even and CP-odd neutral scalar fields, respectively, 
each containing a mixture of Higgs and sneutrino fields. The matrices Z/v, Zr, Za are those 
that diagonalize the neutralino, CP-even, and CP-odd Higgs boson mass matrices and 
are defined in ( |2.22| ) and ( |A-8 , A-1C| ), respectively (for analytic expressions for Zr,Za see 
Eq. (3.14) and (3.25) of [35]). Individually, the neutral fermion - neutral scalar diagrams 
are large, however, if there were no splitting between the mass of CP-even and 



m 

CP-odd neutral scalar eigenstates there would be an exact cancellation between the two 
diagrams. Notice also that the whole contribution is multiplied by a neutralino mass 
which is generically of the order of the electroweak scale. It is rather instructive to simplify 

and m 2 as, 

s>2 A s>2 



Eq. ( |3.5| ) by expanding around 



J Npp 



EE 



--1 i- 



' 4(47r) J 



J V ip 



l 2 



-Z 



c w 



Nlr 



-Z 



S W 



N2r 



Ami, 



m 



2 \2 



hi 



9 



(3.6) 



m 



v+i 



m~ . is the CP even - CP odd sneutrino square mass difference. Its 



where Am?j 

analytical form can be derived from Eqs. (3.14) and (3.25) of [35] to be 

Bf tan 2 



Am?- 



Ml 



Mf 



(3.7) 



where Ma is the CP-odd Higgs mass and Mj the soft breaking slepton masses which are 
diagonal in our basis, see Ref. [35]. A similar expression has been derived in Ref. [18]. Z v 
and Zn are defined in ( 2.18[) and ( |2.22 ), The contribution ( |3.6| ) is driven by the lepton 
number violating terms in the soft supersymmetry breaking sector, B{ and the whole 
expression for the neutral scalar contribution collapses approximately to 



M 



Bf tan 2 P 
{m 2 -M^ 



(3.8) 



where M is the sneutrino or Higgs and m K o the neutralino masses in the loop, respectively. 
The importance of this contribution has already been pointed out in Refs. [15, 19]. The 
mass insertion approximation diagram reads as 



h° 2 



h° 2 



B/W° 



B/W° 



where the 'blobs' indicate insertions of Bs, scalar and gaugino masses. The neutral scalar- 
fermion contribution is thus a ) suppressed from the CP-even-CP odd sneutrino mass square 
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difference, i.e, the lepton number violating soft SUSY breaking parameter Bi, 6) is enhanced 
by tan 2 0, and finally c) suppressed by three powers of SUSY breaking masses. 

The approximate formula ( |3.6| ) does not in general capture the full neutral fermion- 
scalar correction. There are other corrections of the same order of magnitude, including the 
Higgs bosons in s = 1, 2 states. This expansion is more complicated than ( |3.G| ) and is given 
explicitly in section ||, Eq. ( |4.10| ) , where we discuss our numerical results and compare with 
approximate formulae of this chapter. 

Ignoring other than the above diagrams possible cancellations, [Bi tan (3\ must be 
smaller than the 0.1% of the sneutrino mass squared, M 2 , in order to have m u < 1 eV. On 
the other hand, numerically, if the 'solar' neutrino mass difference were to be generated 
by this diagram, then Bi ~ C(l)GeV 2 . Because B { is in principle not constrained from 
above by other means, we conclude that this diagram dominates the whole contribution 
especially when the trilinear couplings, A, A', are negligible. 



3.2.2 Charged fermion - charged scalar contribution 

This contribution reads diagrammatically as: 

H+ 




Using the generic formula the self energy, (C-14), the Feynman rules from Appendix |B.2| , 
and also by applying Eq. ( |3.4| ), we find that the full diagonal contribution for p = q = {5, 6}, 
is : 



J Npp 



EE E E 

s=l r=l i,j,k,l=l a, /3=0 



TO, 



(4vr) 2 



^alk Z H{2+a)s Zl +(2+k)r Z 'N(4+l)p ) x 



Z 



- Z H(2+i)sZ-l r Z N (A+i)p ~ ^/3ij^H(5+j)s^-(2+P)r^N(4+i)p 



B ( 



to 2 o,toL 



, TO 



(3.9) 



where Zjy, Z + ,Z-, Zh are rotation matrices in the neutral fermion, charged fermion, 
and charged scalar sectors, and defined in (2.18), ( 2.31| ), and (A-ll), respectively. It is 
important to notice that following ( p. 31 ) we obtain, Z + ^+k)r — %lkr, with r > 2 and 
hence the contribution ( |3.9| ) is proportional to the mass of a light fermion, m K - . In 
addition, since Z N ^ + ^ p ~ Z u ip, ( P-18| ) shows that the contribution ( |3.9| ) contains the 
rotation mixing matrix Z u , which has been presented analytically in ( |2.22| ). In order to 
analyze the dominant pieces from the charged scalar - fermion contribution, it is instructive 
to consider two cases : A ijfc = and \- k / 0. 

In the case where the trilinear superpotential couplings are absent the charged lepton 
loop has a small contribution to the massless neutrino eigenstates. From the discussion 
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above and we obtain at the limit of small lepton masses (compared to the SUSY 

breaking ones), 



J Npp 



TO, 



^2 7^2 ( X 0lk Z H2s Z lkr Z ulp) X 
s=l i,j,k,l=l r>2 ^ > 



S W 



B {0,m 2 H+ ,0) , (3.10) 



where Ajoj = — Aoij is the lepton Yukawa coupling obtained from Eq. ( |2.27| ), We can 
analyse further equation ( |3.10| ) by Taylor expansion with respect to m' 1 + (commonly named 

^ s 

"Mass Insertion Approximation", or MIA, see, for example, review in Ref. [43]) and using 
( |2~28jg3l|) and ( |A^JaT2|) , 



Z-lr 



Z *H2s Z H(2+i)s m2 H ^ 
Z*H2 S 2 H(5+i)s m2 H -\ 



d 



c z, 



M c 



lr 



JM_ 

M c 



M 



H+2,2+i 



~ Hi<mi tan P , 



M 2 H+ 2, 5+ i ^BitBnP 



(3.11) 

(3.12) 
(3.13) 



where mi is the lepton mass and Mc a generic gaugino mass. Hence, the neutrino couples to 
either the right handed component of the electron, the W~ or the Higgsino with couplings 
proportional to /ij. All the above can be diagrammatically depicted with mass insertions 
as : 



hi 




R 



' M c M 



H+ 



(3.14) 



mi 



Hi mmi 



M c Ml 



tan P 



H + 



(3.15) 



where M^ + is a generic charged Higgs mass. Obviously, both the fermion and the scalar 
propagator are suppressed by lepton number violating couplings. This contribution is then 
compared with the one previously considered with neutral particles in the loop. Indeed, 
in order to account for the atmospheric neutrino mass scale it must be /ij << y/Bl ~ 
0(1 GeV) and hence the charged particle contribution is always smaller than the neutral 
one. Our finding here is in general agreement with the discussion in Refs. [19,27]. Finally, 
notice the Goldstone contribution vanishes since this always conserves lepton number. 

If the trilinear superpotential lepton number violating coupling A^ is turned on, then 
a lepton - slepton loop contribution is generated. In contrast with the pure bilinear case, 
the trilinear contribution may dominate depending on the magnitude of A. In this case the 
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full contribution in (^8|) results in 



yD 
^Npp 



EE E 

s=l r=l i,j,k,l,m,n=l 



m 



(4tt) 5 



\nj Z H(5+j)s Z -(2+n)r Z N(4+l)p 



+ 5 



m 



(3.16) 



Again, making use of (|2T22L ^3ll) we see that the contribution is proportional to the light 
lepton masses and involves the neutrino mixing matrix. We can go a little bit further and 
perform MIA expansion of fl3.16|) as we did before. The contribution then reads, 



J Npp 



mi 



3 

^ ! ITT) 2 

i,j,k,l,m,n=l 



{M 2 



(M 2 



(M 2 H+ )2+m 
~ H +)5+j (M 2 H+ ) 5+j 



■In 



(3.17) 



where mi q is a light charged lepton mass, {M. 2 H+ ) is the charged scalar mass matrix in 



the interaction basis and is given by ( A-12 ). In our notation (A4 H+ )s + j = (Ai H+ )5 + j^ + j, 
and so on. In the denominator and logarithm of ( [3.17 ) one has the difference of diagonal 
elements mm and jj of LL and RR slepton mass matrices, respectively. The approximation 
(3.17) is proportional to the mixing matrix elements {M. 2 H+ )2+ m ,5+j which is nothing other 



than the LR mixing elements of the charged slepton mass matrix ( |A-12| ). These matrix 
elements are (almost) unbounded from experiments when m = j in contrast to the case 
m / j. 

This contribution has been discussed largely in the literature, see for instance [18-21, 
24,27,29]. It is instructive to draw the mass insertion approximation diagram corresponding 
to O : 



'Lm .-< 



'Lp 



-Rk 




'Lp 



Ln 



In the case of dominant \ ikk coupling 



yD 

N 



A m 



2 no tan (3 + Ai 
~W 2 



(3.18) 



with Ai being a trilinear SUSY breaking coupling and M a generic soft SUSY breaking 
mass for a slepton. Comparing ( |3.18| ) with ( p. 14 , 3.15) of the previous case with A — > 
0, we see that the latter is suppressed with at least a factor fii/M. In the case where 
the final two indices are different, A ifc ;,fe 7^ I there is an extra suppression from slepton 
intergenerational mixing and the couplings must be stronger if the lepton in the loop is 
lighter. Our calculation is general enough to allow for these effects too. Furthermore, it is 
obvious from ( 3.17 ) that the r — f-contribution, A i33 , is the dominant one and this coupling 



tends usually to be strongly bounded. 
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3.2.3 Quark - squark contribution 

In general, this contribution originates from up and down quarks and squarks in the loop: 



(la 



II 




Lr 



l Rr 




'Lr 



l Rr 



The up-quark-squark contribution vanishes identically, for the mass eigenstates which 
are massless at tree level. This can be easily seen by applying the master equation ( fTj| ) 
to the corresponding [neutralino-up-quark-up-squark] vertex given explicitly in the Ap- 
pendix p3.3j . 

The case of down quark-squark contribution (the right Feynman diagram above) can 
be divided in two cases depending on the dominance of the trilinear superpotential con- 
tribution : If A' — ► and the only source of lepton number violation is the bilinear term 
then the contribution vanishes. Note that this does not necessarily disagree with the find- 
ings of Refs. [12, 13] where apparently this contribution is claimed to be the dominant 
one. Recall that we are working in the basis of [35] where the sneutrino vevs are zero and 
thus we cannot directly compare, at least graph by graph with this work. In the case of 
Ref. [12] for example, the bilinear term, (iiLiH2 is rotated away. This rotation generates 
new, non-negligible superpotential trilinear couplings which is the case we are about to 
consider. Hence, if \[j k ^ 0, then the situation changes dramatically. Following ( C-14j ), 



the Feynman rules for the down type quarks of the Appendix [B.3| , we find that the most 
general contribution to the massless neutrinos, p = q = {5,6}, reads as, 



J Npp 



E E 



i,j,k,n,m= 



3m d . 



^ (471-) 
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^jikKkm Z } is Z d(3+m)s Z N(4:+j)pZN(A+n)p 



D ( 2 2 2 \ 



(3.19) 



where the rotation matrix in the down squark sector, Zg, is defined in (A-13), and in 
d2l8| ). It is much more instructive to Taylor expand the full contribution ( |3.19| ) around a 
constant SUSY breaking mass into parameters of the original Lagrangian. In the limit of 
small neutrino and quark masses, this results in 



Z v jp Z un p 
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J Npp / 4 \ 2 

j,n,k,m=l 



3 m d k ^jik^nkm 
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In- 



(Ml 



(3.20) 



where the mass matrix M 2 - is defined in (A-14). Notice that (M. 2 ^}i^ +rn are the elements 
of the LR mixing block of Ai 2 ~ and our notation reads (-M|)i = (.M|)jj. The Feynman 
diagram with quark and squark mass insertions representing ( |3.20| ) is: 
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Some remarks are in order : First the quark-squark contribution is proportional to neutrino 
mixing through the matrix ( p. 22 ), and hence to possible hierarchies between fiis. Second, 



it is proportional to squark flavour mixing. Experimental results for K — K, and B — B 
mass difference set severe constraints in the intergenerational squark mixings in the lepton 
number conserving MSSM [(A^|)i,3-f m must be small for i/m]. Although, our calculation 
is as general as possible and allows for these effects we shall assume (AdVfifi+m = 0, i ^ m in 
our numerical results below. The quark-squark contribution may be dominant for sufficient 
large A' couplings. 

3.2.4 Neutral fermion - Z gauge boson contribution 

The corresponding Feynman diagram is : 




Due to the approximate unitarity of the neutrino sub-block of Zn, the contribution of this 
diagram is suppressed either by the lightness of the particle in the loop or by the value of 
the coupling. However, as we obtain from Eq. ( C-15| ) and Appendix ^ this contribution 



is gauge dependent. The dependence again cancels the neutral fermion-scalar contribution 
in ( |3.5[) with the Goldstone boson (s = 1) in the loop. Although we prove this cancellation 
numerically, it can be also shown analytically. 



3.2.5 Charged fermion - W gauge boson contribution 

The Feynman diagram for this contribution is: 




Following (|3.4| ) and the Feynman rules of Appendix B^, when the external legs are purely 
neutrino interaction eigenstates (p, r = 5 or 6), there is no k°-k + vertex. Hence, the 
contribution of this diagram vanishes identically. 
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3.2.6 Summary of the one-loop radiative corrections to massless neutrinos 

The total one loop contribution to massless neutrino masses is given by the sum of the 
neutral scalar loop in (^), the charged scalar loop in ( |3.9[ ), and the squark loop in ( 3.1SQ . 
The gauge boson contributions are negligible. If the trilinear superpotential couplings are 
tiny then the dominant contribution arises from the neutral scalar fermion loop and is 
proportional to CP-even - CP-odd sneutrino mixing [see Eq. ( |3.6| )], If trilinear couplings 
are not small, then depending upon their nature A or A' dominate through lepton— slepton 
[see Eq. (3.16f)] and quark — squark [see Eq. (3.20)] diagrams. 



3.3 Comparison with Literature 

Our work improves on other work which can be found in the literature as no assumptions 
or approximations need to be made. Calculations can be performed in the most general 
supersymmetric model with minimal particle content, without any assumption that matri- 
ces are flavour diagonal, or that any complex phases are set to zero. We have not neglected 
any terms or phases in the neutral scalar sector [35], a basis was chosen in which to perform 
the calculation that had a decoupled CP-odd and CP-even sector and two real vevs. In 
choosing this basis, it is clear that the lepton Yukawa matrix is not, in general, diagonal 
and the lepton mixing matrix does not come purely from the neutrino sector. This is in 
contrast to previous work where, in whatever basis the calculation is performed, the lepton 
Yukawa is chosen to be diagonal. In [14] assumptions are made in the soft sector, such as 
intergenerational mixing being zero, which allows a basis to be chosen where the Yukawa 
matrices are diagonal. Similarly, in [18, 19] there is the assumption of CP conservation in 
the neutral scalar sector. 

Many diagrams are suggested in the literature as being important in generating a 
correct solar mass difference. Under the assumption that the solar mass difference comes 
solely from loop corrections to eigenstates which are massless at tree level, in a general 
basis the external legs must consist purely of neutrino and down-type Higgsino interaction 
eigenstates (in the basis with sneutrino vevs rotated to zero the external legs must consist 
purely of interaction state neutrinos). As such, when diagrams are presented with 'mass 
insertions', it is clear that any diagrams with insertions coupling the neutrino to an up-type 
Higgsino or gaugino on the external leg will not contribute to the solar mass difference at 
one loop. In a basis where sneutrino vevs are not zero, the diagrams with an insertion 
mixing between interaction state neutrinos and the down-type Higgsino contribute to the 
radiative correction of massless tree level eigenstates. In the basis where sneutrino vevs 
are zero this contribution is included in the trilinear vertex, A^. 

Many papers [18-21,27] note the contribution of the loops driven by trilinear couplings 
A^ ) and produce expressions, often with flavour mixing suppressed, that agree with the 
expressions given here. 

The contribution to the charged scalar loop from bilinear couplings is also widely 
noted. Whether a contribution is due to bilinear or trilinear couplings is a basis dependent 
statement [20]. We agree with the results in [12, 14, 19], however in our basis the diagrams 
in [12, 14] are accounted for in the trilinear loops. 
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The importance of the neutral scalar loop has also been noted previously. We agree 
with the general result of [16, 17] that a sneutrino mass difference will give rise to a ra- 
diative correction in the neutrino sector and with [19] that this loop can be the dominant 
contribution. The neutral scalar contribution is included in the analysis presented in [13] , 
but is not discussed in [14]. 

The role of tadpole corrections is stressed in [13] . If we assume the solar mass difference 
arises from the loop corrected 'massless' neutrinos, we can see that the tadpoles do not play 
a role the determining its magnitude. In the interaction picture, there is no ^ Q -^ Q -Higgs 
vertex, so the tadpole contributions vanish. Of course, the tadpoles will affect the other 
heavy neutral fermions. 

A certain class of two-loop diagrams and resulting effects on bounds for lepton number 
violating couplings have been considered [25] 



4. Numerical Results 



In this section we present our numerical results for the neutrino masses. As we have 
already explained, in our most general analysis we use the MNS matrix defined by neutrino 
oscillations as an input. Of course this matrix is not accurately known, but its general 
'picture' has been emerging during the last five or so years with angles and the 3a allowed 
ranges of the neutrino oscillation parameters from a combined, global data, analysis [49], 
reading, 

sin 2 0i2 = 0.24 - 0.40 , sin 2 9 23 = 0.34 - 0.68 , sin 2 13 < 0.046 , (4.1) 
Am^ = (7.1 - 8.9) x 10~ 5 eV 2 , jAm^l = (1.4 - 3.3) x 10~ 3 eV 2 . (4.2) 

In our analysis we fix the neutrino mixing angles to reproduce the tri-bimaximal mixing 
scenario of Ref. [50] , 

sin 2 9 l2 = - , sin 2 6 23 = - , sin 2 13 = , (4.3) 



in agreement with (4.1); the resulting predictions for neutrino mass squared differences are 



then compared with ( |4.2| ), to see whether the values chosen for the input parameters give 
results in agreement with current experimental limits. At present, there is no experimental 
evidence for CP-violation in the leptonic sector; as such, although our analysis is general 
enough to accommodate these effects, in what follows, we shall assume that they are 
negligible. 

In addition to the experimental inputs for the quark and lepton fermion masses and 
mixings, soft supersymmetry breaking masses and couplings must also be initialised. We 
follow the benchmark SPSla [51] where 

M = 100 GeV , M 1/2 = 250 GeV, A = -100 GeV , tan/3 = 10, fi > , 

(4.4) 
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and read the low energy SUSY breaking and superpotential parameters at low energies 
using the code of Ref. [52]. The input parameters of primary interest are those which 
violate lepton number. In the basis of [35], these are, 

Mi j Bi , Xijk , \j k , hijk , h i j k , (4-5) 

where the last two, h and h' are the trilinear lepton number violating parameters in the 
supersymmetry breaking part of the Lagrangian. Apart from these latter parameters, which 
concern trilinear couplings of scalar particles, all others can be used to set the atmospheric 
neutrino mass 2 difference or the solar mass 2 difference. There are two main cases: 

• Tree level dominance : the atmospheric mass 2 difference originates from tree level 
contributions to neutrino masses (Fig. ||). 



1 — loop corrections 




Figure 1: Neutrino mass scales : tree level dominance 



Loop level dominance : The atmospheric mass 2 difference originates from one-loop 
contributions to neutrino masses (Fig. 0). 



1 — loop corrections 




Figure 2: Neutrino mass scales : loop-level dominance 



In either case, the solar mass 2 difference originates from loop effects from the lepton number 
violating parameters in (4.5). 

The correct neutrino mass hierarchy can be always generated by the proper choice of 
just two of the lepton number violating parameters from the list of ( f4.5|) - one of which sets 
the scale of the atmospheric mass 2 difference, the second setting the solar mass 2 difference. 
Of course, in the most general case all parameters can contribute. 
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After choosing the lepton number violating (LNV) parameters, the method described 



in section 2.4 is then employed to determine the charged lepton Yukawa matrix. In general 
it needs to be non-diagonal, in order to reproduce correct masses of both neutral and 
charged leptons and the Umns mixing matrix. The non-diagonal Yukawa matrix (thus 
also non-diagonal charged lepton mass matrix), may easily give rise to effects which are 
already subject to strong experimental bounds; tree level lepton flavour processes, such as 
\i — ► e7 or fj, — > eee, are not suppressed and loop corrections to the electron decays will 
have contributions proportional to the tau mass. To avoid such problems, the specific cases 
considered in the next sections are those for which the large mixing in the lepton sector, 
as seen in the MNS matrix, has its origin purely in the neutral sector, and the charged 
lepton Yukawa couplings remain flavour-diagonal. The formalism we have described thus 
far allows the correct masses and mixing of charged leptons to be initialised. However, 
this will lead, in general, to an off-diagonal lepton Yukawa matrix. These, less natural, 
initial parameters are not necessarily ruled out and within the framework set out above, 
it is entirely possible to perform the calculations as described. However, we now prefer to 
consider a set of parameters for which we do not rely on cancellations in the charged lepton 
sector to make the model phenomenologically viable. The simplest way in which this can 
be achieved, is to find LNV parameters for which lepton Yukawas are diagonal. 



From Eq. ( 2.37| ), for the case where the lepton Yukawa is diagonal and therefore Z\ is 



the unit matrix, we see that 



Zu = U{ INS , (4.6) 



up to higher order terms. Using the MNS matrix as an input, it is possible to see which 
ratios of entries in the mass matrix give rise to the correct leptonic mixing, being, 

mf = Z* diag(mi, m 2 , m 3 ) Z\ = m k U M nsh UjMNSkj • ( 4 -7) 



For example, to set the atmospheric scale at tree level, we can see from Eq. ( 2.20 ) that as 
long as the hierarchy of /Zj matches the ratios of any of the rows in the MNS matrix, the 
mass matrix follow the correct pattern to be consistent with the observed mixing matrix. 
If a second mass scale is set up using the pattern dictated by one of the other rows of the 
MNS matrix, the full MNS matrix will be produced upon diagonalisation. 

It is worth noting the nontrivial fact that such an approach, i.e. generating correct 
structure of neutrino masses and mixings, while keeping FCNC processes in charged lepton 
sector suppressed, is at all possible. 

4.1 Tree level dominance scenario 



At tree level, the mass of a neutrino can be set using m parameters Q2.21 ). The top left 



panel in Fig. |3| shows how the value of (Am^ TM ) varies with \x\ only, setting ^2,3 = 0. The 
grey (or red in color) line is the result given by diagonalising the full neutralino matrix in 
( 2. 14]) and the dark (blue in color) line is given by ( 2.21| ). They agree perfectly in Fig. ||a 



and thus only one line is shown. The shaded band shows the current 3<r limits. From 
Eq. ( |2.21 ) it can be seen, however, that it is |//i| 2 + 1^1 2 + 1^3 1 2 which sets the mass of the 
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All LNV=0 except ^ All LNV=0 except |X| fixed, B 1 varied 




Figure 3: Predictions for atmospheric and solar neutrino mass 2 differences (Am^ TM ) and (Amg OL ) 
for the tree level dominance scenario vs. variations of Lepton Number Violating (LNV) parameters 
as displayed in figure titles. The 3<r gray (green) band consistent with experiment is displayed for 
comparison, as well as the full (gray or red curve) and approximate (dark or blue curve) results, a) 
Only /i! is varied. For all other figures, ^ is fixed as in hierarchy (A) as explained in the text and 
b) B\, or c) A133, or d) A' 133 , is varied respectively. 



tree level neutrino and as such it is straightforward to set any hierarchy between the \i{ and 
still maintain the same value for the atmospheric difference. To correctly reproduce the 
MNS matrix, we choose as an input a very simple hierarchy between the \ii parameters, 

Hierarchy (A) : m = = . (4.8) 

The scale of all three \i{ is set such that a tree level neutrino of the correct mass is generated 
which result in the observed atmospheric mass difference, being, 

Hi = 1.47 MeV , fx 2 = V2x 1.47 MeV , ^ 3 = V3 x 1.47 MeV . (4.9) 

At tree level, this choice of hierarchy gives rise to the MNS matrix, up to the SU(2) rotation 
described earlier, being driven solely by the neutral sector; the charged lepton mass matrix 
is diagonal, and as such we have chosen a set of parameters within this basis which avoids 
the possible phenomenological problems. 

A further, single lepton number violating parameter can then be chosen to set the 
scale of the solar mass 2 difference. The question of the arbitrariness of the tree level 
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neutrino basis is complicated by the requirement that once the loop corrected mass matrix 
is diagonalised, Z\ being the unit matrix is consistent with the experimentally observed 
MNS matrix. As only one further lepton number violating coupling is initiated, the ratios in 
which the loop effects are distributed in the loop corrected mass matrix are approximately 
determined by the tree level mixing matrix. As such we can determine an approximate 
expression for the extra contribution to the full rotation matrix from rediagonalising the 
loop effects. The further condition that the full rotation must reproduce the MNS matrix 
allows us to fix the tree level basis. 

The three further Figs. |3](b,c,d) show the range of possible parameters in this scenario. 
In each of these plots, the set of \ii are given the values fl4.9p and another, single lepton 
number violating coupling is varied. In each case, the gray (or red) line shows the full result 
and the dark (blue) line is the result predicted by the approximate solutions. The fact that 
A133 and A' 133 give the correct solar mass difference over a similar range of parameters is a 
numerical coincidence. For this example, the factors from the different fermion masses in 
the quark loop, colour counting and scalar mixing cancel each other. 

The contribution of the neutral scalar loop discussed in section ^, results from cancel- 
lations between the CP-even and the CP-odd diagrams and may includes contributions of 
approximately the same order. As such, the approximation presented earlier in the text, 
Eq. (|3.6| ), does not agree well with the full result. The discrepancy between the full result 
and the approximate result reflects the fact that various contributions arise from different 
places in the full calculation (e.g. the effect on the mixing matrices, the effect on the 
sneutrino masses). The approximate result plotted in Fig. ||]b is given by 
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(4.10) 



The approximate result for the charged scalar loop, given by Eq. ( p. 17 ) agrees well with the 
full result (Fig. |3|c). However, as A133 = — A313, there are other diagrams which contribute 
to the full result which are not included in the approximation. The approximate expression 
captures the important effect. The agreement between the full result and the approximate 
result, given by Eq. (|3.20|) when varying A' 133 (see Fig f|d) is very good, as the diagram 
highlighted in the text is the only diagram which contributes. 
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All LNV=0 except X i33 All LNV=0 except ^33 fixed, V h1 varied 




Figure 4: Same as in Fig. 3 but for the loop-level dominance scenario. All LNV parameters arc zero 
apart from a) Ai33 that is varied in hierarchy (D). For all other figures, A^33 is fixed to a value (see 
the text) consistent with the atmospheric mass 2 difference and b) only X' in is varied in hierarchy 
(B) or c) only X' i22 in hierarchy (B) or d) only A^ 33 in hierarchy (B) in order to accommodate the 
solar mass 2 difference. 

4.2 Loop level dominance scenario 

It is possible for both the solar and the atmospheric scales to be set by loop corrections. 
This happens if the bilinear parameters fii are small enough. In this section we analyse 
this case setting strictly \x\ = ^2 = M3 = 0> so that the one-loop corrections to the full 3x3 
neutrino mass matrix are finite. Otherwise a more involved renormalisation scheme has to 
be implemented. 

Again, we would like to set the Lagrangian parameters such that one can generate the 
correct structure of the MNS matrix while keeping the charged Yukawa couplings flavour- 
diagonal. This can be achieved if the neutrino mass hierarchy is governed by the trilinear 
A and A' couplings. For the diagrams dominated by trilinear couplings the flavour of the 
external legs of the loop can be "swapped independently" of the flavour of the particles in 
the loop, just changing the appropriate indices of the A, A' matrices in the loop vertices. 
Setting the A and A' entries which control the couplings of the external legs in certain 
hierarchies, one can ensure that also the ratios of the various entries in the one loop 
corrected neutrino mass matrix are such that they give rise to the correct Umns rotation 
matrix. 
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All LNV=0 except V i33 
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Figure 5: Same as in Fig. 3 but for the loop-level dominance scenario. All LNV parameters are zero 
apart from a) A^ 33 that is varied in hierarchy (B). For all other figures, A^ 33 is fixed to a value (see 
the text) consistent with the atmospheric mass 2 difference and b) only X' in is varied in hierarchy 
(C) or c) only A- 22 in hierarchy (C) or d) only A i33 in hierarchy (D) in order to accommodate the 
solar mass 2 difference. 



The possible hierarchies are given by the rows of the MNS matrix and are, with a 



generic coupling A^A as follows 8 , 
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(4.11) 

(4.12) 
(4.13) 



Due to the antisymmetry of the first two indices of A it can only be chosen to follow 
hierarchy (D) described above. 

As the nature of the loop corrections due to B{ means that the external legs cannot 
be swapped without affecting the flavour structure inside the loop, it is difficult to fix a 
hierarchy of Bi in the Lagrangian which will automatically give rise to the correct ratios 
in the one loop corrected mass matrix. 



8 Couplings with A^A, with j ^ k have only negligible contributions to neutrino masses and excluded 
from our hierarchy list. 
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We consider first, the case where the atmospheric mass 2 difference is set by Aj 33 in 
hierarchy (D). The range of values for which the correct atmospheric mass difference is 
given is plotted in Fig. [|a. Note that although we plot on the x-axis A133, the coupling 
A233 is also varying to keep the hierarchy (D) fixed. The fact that both A133 and A233 
contribute is the reason the value of the coupling is only a little greater than the value 
of A133 which correctly reproduces the solar mass difference in the tree-level dominated 
scenario. 

The further three panels [Fig. ||(b-d)] have a fixed set of Aj 33 in hierarchy (D) giving 
the atmospheric difference. Being, 

fi 7 x lfP 5 

A133 = 6.7 x 1(T 5 , A 233 = -—^7| . A 333 = . (4.14) 

A333 = due to the antisymmetry between the first two indices, fitting hierarchy (D). In 
addition to this, Fig. ||]b varies A' m 2ll 311 in hierarchy (B) giving the solar mass 2 difference. 
Again, we plot the solar mass 2 difference against the value of A' m , however A^j 3 u are 
being varied at the same time. The two remaining panels take in turn different sets of three 
A' couplings, \' i22 and A^ 33 . The final two indices determine which particle is produced in 
the loop. As such, with a lighter particle in the loop, the couplings must be greater to 
compensate. We see that, in moving from one panel to the next Fig. |](b^ d), to produce 
the same mass difference, a smaller value of the coupling is required with a heavier particle 
in the loop. With the down quark in the loop (Fig. |p) the value needed for the coupling 
may result in large contributions to the neutrinoless double beta decay rate as it is already 
approaching the excluded regime [53]. 

Finally, Fig. |H|a, we show how the atmospheric mass 2 difference can be set by the three 
A^ 33 couplings in hierarchy (B), plotting the result for the atmospheric mass 2 difference 
against A' 133 . Next, we set A^ 33 to take the following values 

A' 133 = 3.25 x 10~ 5 , A2 3 3 = V2 x 3.25 x 1(T 5 , A 333 = V3 x 3.25 x 1(T 5 . 

(4.15) 

The remaining three plots, Fig. |||(b,c,d), show the change in solar mass 2 difference, as sets 
of either A' in i22 in hierarchy (C) or the set Aj 33 in hierarchy (D) are varied. 

5. Conclusions 

An increasingly accurate picture of the neutrino sector, with masses much smaller than 
the charged leptons and a distinctive mixing matrix in the W-vevtex, is being discerned 
by current experiments. We note that there are three preferred Zn symmetries in the 
supersymmetric extension of the Standard Model with minimal particle content. Imposing 
a Z2 symmetry results in the widely studied R-parity conserving MSSM, however another 
preferred symmetry, Z3, gives rise to a Lagrangian which explicitly violates lepton number. 
These interactions lead to neutrino masses, both through a 'see-saw' type suppression at 
tree level and through radiative corrections. We have considered the most general scenario 
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in this model; no assumptions have been made concerning CP-violation or intergenerational 
mixing, for example. 

We present, in Appendix |A|, the tree level mass matrices of the model and, in Ap- 
pendix |B], the full set of Feynman rules for the neutral fermion interactions. The calculation 
has been performed using two-component Weyl fermion notation; Appendix |C] contains a 
derivation of the propagators in this notation, included primarily for pedagogical reason, 
and generic self-energy diagrams for scalar and gauge boson corrections to fermions. 

In the basis set out in our previous work [35] , we find that a non-zero neutrino mass will 
arise at tree level unless all fii are zero and analyse in detail, the further contributions to 
masses that come from loop corrections. We show that the magnitude of the contributions 
due to neutral fermion loops, examined in section 3.2.1 are determined by the size of 
the bilinear supersymmetry breaking parameter, Bi\ that loops with charged fermions, 
described in section [3.2.2| , have a contribution due to trilinear lepton number violating 
couplings in the superpotential, Ayjt; and that quark loops, section 3.2.3 , are determined 
by the trilinear lepton number violating coupling A^- fe . Each of these contributions can be 
dominant. In sections 3.2.4 and 3.2.5, we consider the gauge loops and why they do not 
give large contributions to neutrinos which are massless at tree level. We derive expressions 
for the full calculation, which form the basis of our numerical analysis. We also present 
approximate expressions in each section, which are simple, compact formulas encoding the 
important information pertaining to each diagram. In our presentation of the results, as 
seen in Figs. |3[|],|| these simple expressions are shown to be in good agreement with the 
full result. 

The lepton sector in the 1/i-MSSM is much more involved than the lepton number 
conserving MSSM. Mixing between leptons, gauginos and higgsinos ensures the question 
of initialising Lagrangian parameters must be carefully considered. A framework has been 
constructed in section ^ to correctly reproduce the charged lepton masses and MNS matrix 
for any set of lepton number violating couplings. 

In constructing the framework in which to perform the calculation it is clear that there 
will be, in general, large intergenerational mixing in the lepton Yukawa matrix, this allows 
the possibility of unsuppressed tree level flavour violating processes, already subjected to 
strong bounds. To circumvent this problem we considered sets of Lagrangian parameters 
for which the MNS matrix has its origin solely in the neutral sector, the lepton Yukawa 
matrix being diagonal. The three rows of the MNS matrix correspond to three sets of ratios 
between entries in the loop corrected mass matrix which will give the correct MNS angles. 
It is possible to set these ratios by setting hierarchies in the couplings between generations. 
With the condition that it must be possible to change the flavour of the external legs of 
the diagram without affecting the flavour structure of the loop, there is some freedom in 
choosing which group of Lagrangian parameters we set in each hierarchy. 

Lepton number conserving parameters were fixed to be the SPSla benchmark point, 
and we have investigated the effect of varying the lepton number violating couplings, as 
seen in Figs. HJ4],[|. We have shown that values for lepton number violating couplings exist, 
which give the correct atmospheric and solar mass 2 difference, charged lepton masses and 
mixing, which are not already excluded by existing studies of low energy bounds. There 
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are two distinct scenarios that achieve this: the tree level dominance scenario, in which 
the atmospheric scale is set at tree level and the solar scale set by radiative effects, and 
another, the loop level dominance scenario, in which both the atmospheric and solar scales 
are set by radiative corrections. 

In the tree level dominance scenario, we choose the /ij parameters to be of the order of 
1 MeV, such that the correct result for the atmospheric mass 2 difference is obtained. They 
are chosen to obey a certain hierarchy, which ensures the mixing matrix is consistent with 
observed MNS. 

In addition to this, a further, single lepton number violating coupling can set the 
scale of the solar mass 2 difference by determining the contribution of the appropriate loop 
diagram. It is possible to generate loop diagrams of the appropriate scale, by including 
either a non-zero A, A' or B coupling. We find that the correct solar scale can then be set 
by any of 

B x ~ 0.21 GeV 2 ~ [300 Aii 
A133 ~ 3.4 x 10~ 5 ~ y e , 

A' 133 ~ 3.2 x 1(T 5 ~ 0.1y d , (5.1) 

where y e ^ is the Yukawa coupling of either the electron or the down quark, presented here 
merely for the sake of comparison. 

In the second case, the correct masses and mixing for both charged and neutral fermions 
can be achieved without a massive neutrino at tree level. The solar and atmospheric mass 2 
differences both arise from radiative corrections at one loop, using loop contributions whose 
value is determined by sets of A or A' couplings in given hierarchies, such that the observed 
MNS is generated. Firstly, we show that we can set the atmospheric scale with a set of A 
couplings of the order of the electron Yukawa coupling, then find the solar scale is correctly 
set by A' couplings of the order of the down quark Yukawa coupling. 

Alternatively, the atmospheric scale can be set by A' couplings, 

Ai 33 = = = 3-25 x 10- 5 ~ 0.1y d , (5.2) 
and the solar mass 2 difference can be generated by another set of A' couplings, 

A' 122 = = -^f ~ 6.5 x ID" 4 ~ 2y d , (5.3) 

or a set of A couplings of the order of the electron Yukawa. 

We include some comments on how this work compares with previous work in the 
literature. We highlight where our results agree with statements made in the literature 
and comment on results presented in different bases. 

The study of neutrino masses will provide the basis for further work concerning lepton 
number violating phenomena. The ranges of values for lepton number violating parameters 
required to produce the correct masses and mixing, will be reflected in processes such as 
tree level lepton flavour violating decays and will have repercussions concerning rare events 



- 32 - 



such as neutrinoless double beta decay. This will make a valuable link between collider 
experiments and upcoming neutrino experiments. In this paper, we have made a framework 
for these investigations. 
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Appendix 



A. The Lagrangian and the mass matrices of the ^,-MSSM 

We strictly follow the notation of Refs. [23,35]. The L'-MSSM superpotential is given 
by Eq. (2.1). The main discussion of this paper is confined in the lepton sector, but for 
completeness we define here also the mass basis in the quark sector. To this end, we rotate 
the four quark superfields to a basis where both \' 0i j and (Yu)ij are diagonal 

D' L — Z* dL D L , D' — Z dR D , 

U' L ^Z UL U L , u'^z: R U. (A-l) 

By absorbing a rotation matrix into the Lagrangian parameters, one can write down the 
superpotential Q2.1|) as 



+ Oru)ij Q} X H 2 2 U* - (Yu) kj K ki Qf x H 2 1 Uj , (A-2) 

where X' j k and {Yjj)ij are diagonal matrices and C 1 (£ 2 ) is the neutrino (electron) com- 
ponent of the SU(2) doublet. The charged current part of the Lagrangian, 

C w = -^—ULiKij^W+dLj + H.c , (A-3) 
diagrammatically reads (in Weyl spinor notation) as, 



V2sw 5 



(A-4) 



with Kij = Z* 

ki^d L kj being the CKM matrix. We rotate all fields in the basis where 
sneutrino vevs are zero. Then the soft supersymmetry breaking terms are, 

£ssb = - (M?)^ v* La v L/3 - (M?)^ e* La e Lp - m 2 H2 hfh° 2 - m 2 H2 h+*h+ 

~ (h^ijuiih^ugj + ih u ) kj K kl d y Li h^u^ - h af3k D La e Lf3 e Rk - h' ajk u La d Lj d Rk 
+h' ak jK* k e La u Li d* R j + B a i> La h2 — B a e La h\ 



- 34 - 



+ -MiBB + M 2 W + W~ + ^M 2 W°W° + ^M 3 G R G R + H.c. 



(A-5) 



where B a is the four-component bilinear term B a = (Bo,Bi) and h, h' are trilinear soft 
breaking couplings. In the basis of [35], in addition to vanishing sneutrino vevs, one obtains 



also diagonal sneutrino soft breaking mass terms, ( M~ 



M~ is given in Eq. (2.25) 



of [35]. Notice also that the soft breaking mass which corresponds to the mixing of the Higgs 
with the slepton is just the term Bi tan (3. In this basis, we shall now present the spectrum 
of the model. The neutral Higgs sector and approximate formulae has been displayed in 
Ref. [35]; we repeat only the mass matrices here for completeness and definition. 

A.l Mass Terms for CP-even Neutral Scalars 

After electroweak symmetry breaking, sneutrinos, {vii}, mix with Higgs bosons, {Plo = 
hi,h 2 }, resulting in CP-even and CP-odd scalars (recall that CP-symmetry is preserved 
at tree level even in the most general R-parity violating MSSM [35]). The corresponding 
CP-even neutral scalar mass terms are 
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The rotation matrix is then given by 

Z\M\Z R = diag[mlo,m 2 H o,(ml + )i] . i = 1, ...3 . 
An approximate formula for the matrix Zr is given in Eq. (3.9) of Ref. [35]. 

A.2 Mass Terms for CP-odd Neutral Scalars 

Mass Terms for CP-odd Neutral Scalars can be read from 



Im h 2 , 



C D 



(im/i 2 , ImP L0 Im i>Li) z a Z A<M.\Z ' a Z\ ImP L0 

lmi> Lj 



where 
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Bi 



Bi tan (3 



M 2 Sij j 



(A-6) 

(A-7) 
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(A-9) 



- 35 - 



and £ is the gauge-fixing parameter. The rotation matrix Za is defined through, 

Z T A M 2 K Z A = diag[m 2 G o,rn 2 A o,(rnl_)i] , % = 1, ...3 . (A-10) 
An approximate formula for the matrix Za is given in Eq. (3.22) of Ref. [35]. 
A. 3 Mass Terms for Charged Scalars 

In t'-MSSM charged Higgs and charged sleptons mix. The mass terms and the rotation 
matrix Zh, can be read from the Lagrangian 

(ht\ 

£D-(/i 2 + * h~ x e Lj e Rk )M' H , 



-Li 



(hp K hj ^Rk)zH^ H M 2 H+ Z H Z^ H 

\ e ra / 

= -(M 2 H+ ) pq H+*H+ , p,q=l, ...8 , 
where the notation is self explanatory and 
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M\ sin /3 cos (3 + (1 - £) sin /? cos /? M| sin 2 /3 + cos 2 /3 + sin 2 (3 



H+ 



Bj 



V 



Bj tan /? 



M 2 h 



Bi tan P 73 A 0m«/4^ 
- M 2 , cos 2 2/? + iA* im A 0im t; 2 ^( - \ ajl fi* a v u + /^d) 

^ ( - KjlWu + h* 0jl v d ) (m| c ) ^ + (M 2 , - M 2 ) cos 2 2/? % + \\ Qml X* mk v 2 ) 

(A-12) 

A. 4 Mass terms for down-type squarks 

Mass terms and rotation matrices for down-type squarks arise from the Lagrangian part 



-C^(dfJ^)M 2 A 
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Rl 
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l Rl 
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(A-13) 



where 



M 2 



m 



V2 



(A-14) 



Recall that A' ofcm = Yok dkm are diagonal down-quark Yukawa couplings. 
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A. 5 Mass terms for up-type squarks 

The mass terms for up-type squarks are 



£ = ( fig ) -^u ( ~| J = ( as e% ) ^4^4 



'Life 

2 



^2 ~*2~2 



upg ^p ) 



p,g = 1, ...,6 . 



(A-15) 



where 



mi 



Km l K] ) lk + h^uY^uvl + (i - 9 f)(v 2 u - v 2 )S ik j= (K) jk v u - (Yu)* k v 

73 {hu)u vu-^ {Yu) u ii* v { m h) jm + U y u Y u hmvl ~i(v 2 u - v 2 )S jm 



Recall that (Xu)ij = Yui $ij are diagonal up quark Yukawa couplings. 
A. 6 Mass terms for down quarks 



7T 



£ ^ -—/^^Oij v dd Z Lid Z Rj - -^^Oij v dd Z LidRj 



= -m di d z Li d z Ri - m di d z Li d Ri , i = 1, ...3 . 
A. 7 Mass terms for up quarks 

£ D ( Yu h VuU Rj U Li ~ ^ ( Y uTij V uU V Rj U V Li 

= ~ m uiU Ri u y Li — rn u iU y Ri u y Li , i = 1, ...3 . 
A. 8 Mass terms for neutrino-neutralino 
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where 
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A. 9 Mass terms for charged lepton-chargino 
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where 



B. Feynman Rules for Neutral Fermions in the ^-MSSM 



B.l Neutral Scalar - Neutral Fermion - Neutral Fermion interactions 
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B.2 Charged Scalar - Neutral Fermion - Charged Fermion interactions 
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B.3 Squark - Neutral Fermion - Quark interactions 
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B.4 Fermion - Fermion - Gauge boson interactions 
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C. Weyl spinors and self-energy one-loop corrections 



Throughout this article we make an extensive use of Weyl spinor notation. Generally 
speaking, working with Weyl spinors is advantageous because of two reasons : first they 
appear naturally in a supersymmetric Lagrangian (no extra work is required to make 
the connection with Dirac or Majorana four-component spinors and their corresponding 
Feynman rules) and second, when used in Feynman diagrams, they present transparently 
their structure as for example, the dominance of a particle mass or the appearance of a 
mixing etc. In this Appendix we give a pedagogical introduction to the use of Weyl fermion 
propagators and vertices. We then calculate at the end the generic self energies that appear 
in (3.1, |3.2| ). This Appendix is complementary to the work of Ref. [45]. 

Using path integral technics we want to find the propagator of a massive Weyl fermion 
in Minkowski space. The path integral involving a Weyl fermion £ a (x) and a source field 
J a {x) is 

W[J,J] = N j d[^} d[^} e iSl ^' J ' J] , (C-l) 
where N is a constant and 

i S[£, £, J, J] = i J d 4 x| \ ila»d£ + itaHdJ - m(tt + ff) + J£ + Jf| . (C-2) 

In writing the above action functional we use the conventions of Wess and Bagger [48] with 
the metric being g^ u = (1, —1, —1, —1). It is simpler to work in momentum space and thus 
we Fourier transform (in four dimensions) both fields and sources as 

/•+°° d A v • -■ f +0 ° d^n 

rw = L w? e "" e(p) • eix) = L w? e "" r(p) ■ 
= C w? e "" J " ip) • Ji{x) = C w? e " r " J4(p) • (a3) 

where p ■ x = p^x^. We also make use of the four dimensional definition of the 5-function, 

r+co i4 

5{x-x') = / 77^4 e*H*-*') . (C-4) 



(2vr) 4 

The action functional ( |C-2| ) is conveniently written in a matrix notation 

d 4 P ~ (n^Mn + rtx + x^nj , (c-5) 

where 



fi(p)= I. , X(p)= , M(p)= , (C-6) 
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where a ■ p = -p u and a ■ p = at an d the matrix A4 is hermitian. The path integral 
measure does not change when transforming the Weyl fermions £ by a constant (J) as 

n = n' - at 1 x . (c-7) 

After this transformation is applied in ( p-5 ), sources and fields "get decoupled" 

i J d 4 p^ftt^Mn' - X^M^X) , (C-8) 
where the inverse of M is easily found to be 



AT 1 



p2 _ m 2 



a ■ p m 
m a • p 



(CM)) 



The first integrand term in ( |C-8| ) is exactly the same as the one in square brackets of ( |C-2[ ) . 



Recalling that X = X[J, J] the path integral in (C-l) takes the form 



W[J,J] = Ne iSo ^ exp<j ~ I dtpX^M^X j> = W[0] exp <j ~ / cfpX^M^X 



where So is the free Weyl fermion action functional 
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i^d^ - -m(ti + ft) 



S [U] = 

The propagators for the Weyl fermions can be read from 



(C-10) 
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2 Jp{-p) + U-p) 2 \ J $ (P)\ • (C-12) 

pZ _ m z f pZ _ m z 



Diagrammatically the propagator of a massive Weyl fermion is depicted in Fig. []. The 
convention we adopt here is that arrows run away from dotted indices at a vertex and 
towards undotted indices at a vertex. As it is obvious from (C-12), the kinetic part of the 
propagator [the top one in Fig. ^| is uniquely defined from the height of the indices that 
link this propagator with the vertex. 

The propagator for two different Weyl spinors rj and ip (forming a Dirac spinor in four 
component notation) with action functional, 



i S[£, f , J, J] = i / d 4 x 



+ J V T] + J n ri + Jyip + J^V 
(C-13) 
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Figure 6: Massive Weyl spinor propagators. 
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Figure 7: General vertices involving Weyl spinors used in this article. In case of the vector boson 
vertex the Feynman rule is defined completely from the height that link a propagator with this 
vertex. Vertices on the left arise from Cy = A rst 4>tipsVr + J ^l s t4't' l Ps'rjr an d vertices on the right 
arise from C = B rs 'tj} s a IJj V l _ l r] r + B^fjra^V^tps. 



have the same form as in ( C-12 ) with obvious Lorentz invariant substitutions. Vertices 
with Weyl fermions arise from the superpotential and from the super symmetric gauge 
interactions. They are displayed in Fig. [j]. 

We now proceed in calculating the general self energies of fl3.1|)|3^ ) . The 1PI self energy, 
Ti D , obtains corrections from diagrams which have either gauge particles or scalar particles 
in the loop. The scalar contributions for general scalar-fermion vertices iA qrs and iB prs is 



At v R 



S £( m *g) 




ITT 



iBprsAqrsm^rjr ^U^- D B ° ' m <t>s ' m V"?r ) ' 

(C-14) 



where denotes the physical mass of the mass eigenstate which is composed out of 

the interaction eigenstates ip r and r] r . The corresponding neutrino self energy arising from 
vector boson contributions with generic vertices iCqrU 11 and iD qr a^ is 




(C-15) 



where £ is the gauge fixing parameter, my is the mass of the vector boson and Bq, Co are 
the Passarino-Veltman functions [46] in the notation of Ref . [47] , 



B Q {q .m^tfty) 
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(C-16) 
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(C-17) 

Finally, self energy corrections to the Weyl fermion kinetic terms read as 
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(C-18) 



and 
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-(£+l)B (m 2 K o, my, ml r )-2Bi(m 2 K o, m\, m^ r ) 

- (£- l)£rn v C {0,mlo,mlo,m v ,£m v ,ml r ) 

- (f- 1) (m 2 K o ~m\ r ) C 2 (0, mlo , m 2 K o , m 2 v , im\ , m\ r ) 

(C-19) 

where -Bi , C2 are defined in [47] . 
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